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ABSTRACT 

Results are presented from a long-term astrometry program on PSR B0919+06 using the 
NRAO Very Long Baseline Array. With ten observations (seven epochs) between 1994-2000, we 
measure a proper motion /i Q = 18.35±0.06 mas yr _1 , [i$ = 86.56±0.12 mas yr _1 , and a parallax 
7r = 0.83 ± 0.13 mas (68% confidence intervals). This yields a pulsar distance of 1.21 ± 0.19 kpc, 
making PSR B0919+06 the farthest pulsar for which a trigonometric parallax has been obtained, 
and the implied pulsar transverse speed is 505 ± 80 km a" 1 . Combining the distance estimate 
with interstellar scintillation data spanning 20 years, we infer the existence of a patchy or clumpy 
scattering screen along the line of sight in addition to the distributed electron density predicted 
by models for the Galaxy, and constrain the location of this scattering region to within ~ 250 
parsecs of the Sun. Comparison with the lines of sight towards other pulsars in the same quadrant 
of the Galaxy permits refinement of our knowledge of the local interstellar matter in this direction. 



Subject headings: astrometry — ISM:general — pul 
1. Introduction 

Since the discovery of pulsars, it has been rec- 
ognized that the dispersion and scattering of pul- 
sar signals provide unique information about the 
intervening medium. The distance to most pul- 
sars is estimated using the observed dispersion 
measure (DM) and a model for the Galactic elec- 
tron density distribution (e.g. Taylor & Cordes 
1993, hereafter TC93). Where available, a model- 
independent distance from annual trigonometric 
parallax provides crucial calibration information 
for this model, as well as allowing the absolute 
luminosity of the pulsar to be derived. 

Pulsar astrometry also provides the observa- 
tional evidence required to investigate several 
other questions. Proper motion measurements 
(especially in conjunction with reliable distance 
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estimates) allow verification of pulsar-supernova 
remnant associations (e.g. Kaspi 1998). The 
inferred speeds of pulsars constrain the mini- 
mum asymmetry in supernova core-collapse pro- 
cesses, or other sources of kick velocities. Pulsar 
population statistics, selection-effect biases, and 
planetary-extragalactic reference frame ties have 
also been addressed by astrometry (e.g., Lyne & 
Lorimer 1994; Cordes & Chernoff 1998; Bartel 
ct al. 1996). 

At present, there are only a handful of model- 
independent distances to pulsars. In a recent sum- 
mary, Toscano et al. (1999a) list nine such objects. 
In §2 of this paper, we present results from a long- 
term astrometry program on PSR B0919+06 using 
the NRAO Very Long Baseline Array (VLBA). We 
measure the proper motion and trigonometric par- 
allax for this pulsar, and infer a distance of ~ 1.2 
kpc and a transverse speed of ~ 500 km s , com- 
parable to the mean population speed. 

PSR B0919+06 has a DM of 27.31 pc enr 3 
(Phillips and Wolszczan 1992), yielding a mean 
electron density of 0.023 cm -3 for the interstellar 



medium (ISM) along this line of sight. However, 
interstellar scintillation (ISS) and scattering of the 
pulsar signals can be used to extract much more 
information about the line of sight than simply 
this mean electron density. In §3, we review the 
necessary formalism, and analyze the published 
data spanning ~ 20 years for this pulsar. Our 
analysis combines data from very long baseline in- 
terferometry (VLBI) and ISS, in order to infer the 
distribution of scattering material along the line 
of sight. We find the need for extra scattering ma- 
terial in addition to the TC93 distribution; with 
the assumption that this material is present in the 
form of a screen or clump with variable scattering 
strength, we constrain the scren location to the 
local interstellar medium, within <~ 250 parsecs of 
the Sun. 

This analysis, along with DM and distance mea- 
surements for other pulsars in the third quadrant 
of our Galaxy, allows refinement of our knowl- 
edge of the local electron density distribution. Be- 
sides being a situation of interest in its own right, 
this analysis also serves as a model for the hybrid 
VLBI-ISS technique applied to the local ISM. 

2. Obtaining a VLBA Proper Motion and 
Parallax 

We undertook observations of PSR B0919+06 
using the VLBA, over seven epochs (ten obser- 
vations) from 1994 to 2000. The details of data 
acquisition, processing and calibration have been 
discussed by Fomalont et al. (1999), who pre- 
sented a preliminary result based on the first four 
epochs of this dataset. For completeness, we sum- 
marize the data reduction procedure here. 

2.1. Phase Referencing and In-Beam Cal- 
ibration 

In VLBI observations, the visibility phase has 
to be estimated for each baseline ( "fringe fitting" ) 
in order to correct for variations in clock off- 
sets, station positions and atmospheric propaga- 
tion effects. However, fringe fitting is not possi- 
ble for weak sources, which lack sufficient signal- 
to-noise ratio (SNR) on the short timescales on 
which the visibility phase varies. Self-calibration 
enables corrections of unmodeled phase errors, 
but destroys absolute positional information. As 
such, astrometric observations must use phase- 



referencing, where scans are alternated on the 
(possibly weak) target and a strong nearby cali- 
brator (the "nodding calibrator"). Phase connec- 
tion across scans on the nodding calibrator allows 
the visibility phase to be interpolated between the 
scans on the target, extending the effective coher- 
ence time from minutes to hours. This allows as- 
trometry of weak sources, relative to the (assumed 
fixed) position of the nodding calibrator. This 
scheme requires limits on both the angular sep- 
aration between the target and the nodding cali- 
brator, and the cycle time between them. These 
requircmnts are outlined, for example, by Wrobel 
at al. (2000), who recommend cycle times rang- 
ing from 30 seconds at 43 GHz to 300 seconds at 
1.4 GHz, with an angular separation < 5°. 

The negative spectral index of pulsars can be 
exploited by observations at lower frequencies, 
where they are stronger. However, the achieved 
resolution is proportional to the synthesized beam 
size, and inversely proportional to the SNR. At 
lower frequencies, the beam is larger, and the sys- 
tem temperatures are generally higher, offseting 
some of the gain in pulsar signal. Additionally, 
ionospheric effects become more troublesome at 
lower frequencies. Our chosen compromise was to 
observe at L-band (1.4-1.7 GHz), with a typical 
target-calibrator separation of 5° and a cycle time 
of 5 minutes. Even with this strategy, the residual 
ionospheric errors after phase-referencing can oc- 
casionally be large enough to prevent correct phase 
interpolation and useful astrometry. 

Various techniques have been explored to cali- 
brate out differential ionospheric effects, including 
the use of GPS data (e.g. Ros et al. 2000). For 
PSR B0919+06, we use an in-beam calibrator, as 
discussed in detail in Fomalont et al. (1999): there 
is a faint (~ 10 mJy) source, J0922+0638, 12' from 
the pulsar, within the ^30' primary beam of a 
VLBA antenna at 1.4 GHz. The pulsar data are 
correlated twice, once at the pulsar position and 
once at the position of the in-beam calibrator. We 
phase-reference the data to the nodding calibrator, 
self-calibrate the visibility phases so that they are 
consistent with the visibility function expected for 
the in-beam calibrator, and transfer these correc- 
tions to the pulsar data to carry out the imaging 
and astrometry. 

This procedure has two crucial advantages over 
ordinary phase-referencing. Phase errors from 
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the different ionospheric electron column densi- 
ties along different lines of sight generally increase 
with angular separation. Using an in-beam cali- 
brator not only reduces the target-calibrator an- 
gular separation from ^5° to a few arcminutes, 
but the visibility phase corrections are also derived 
for the actual times when the pulsar was observed, 
and thus avoid interpolation over time. This pro- 
cedure reduces the residual ionospheric errors so 
that sub-milliarcsecond astrometry is achievable. 

2.2. Imaging and Astrometry 

The VLBA data used here were obtained as 
part of two different programs, using different nod- 
ding calibrators and slightly different frequency se- 
tups. The earlier observations (f994-f996) used 
VLA calibrator 0906+0f5 (6° away from the tar- 
get) as a nodding calibrator, and a cycle time of 
4 minutes on target and 2 minutes on the cali- 
brator. The later observations (1998-2000) used 
J0914+0245 (4° away) as the calibrator and short- 
ened the cycle time to 3 minutes on target, 2 min- 
utes on calibrator. We also adopted a strategy 
of snapshot images at two closely spaced epochs 
(3 to 7 days) in order to verify the consistency 
of the astrometry, and slightly adjusted the spe- 
cific observation frequencies (8 independent chan- 
nels of 8 MHz each) within the 1.4-1.7 GHz band 
to avoid some known radio frequency interference. 
During some later epochs, the pulsar was gated 
at the VLBA correlator, boosting the SNR by 
~ Z^ 1 / 2 , typically a factor of 3-4, where / = 
T on / (Ton + Toff) is the gate duty cycle. 

We note that the changes in observing setup 
between epochs do not affect the astrometry, since 
the final pulsar positions are determined relative 
to the (fixed) position of the in-beam calibrator, 
which was assumed to be the same at all epochs. 
The absolute position of the in-beam calibrator, 
determined with respect to J0914+0245, is known 
with a precision of ~10 mas. 

The data were reduced using standard VLB 
phase-referencing procedures (Beasley & Conway 
1995) using AIPS, the Astronomical Image Pro- 
cessing System. This involved amplitude calibra- 
tion using the system temperature at each antenna 
after flagging bad records; a priori phase calibra- 
tion using pulse calibration tones at the VLBA 
when available; fringe-fitting to the nodding cal- 
ibrator; and self-calibration of the in-bcam cal- 



ibrator. The calibration was transferred to the 
pulsar, and both the in-beam source and the pul- 
sar were imaged. During epochs with a relatively 
undisturbed ionosphere, the in-beam source was 
observed to be point-like before self-calibration. 
Self-calibration yields a point source image for the 
in-beam source, whose position is used for relative 
astrometry of the pulsar. Figure 1 shows one pul- 
sar image from each of the seven different epochs, 
mosaiced together to illustrate the proper motion. 
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Fig. 1— Proper motion of PSR B0919+06, de- 
picted as images for a series of epochs. Contours 
are shown at 5,10,20,40c, where a is the RMS 
map noise. There are no pixels in the map at or 
below — 5cr. The period-averaged pulsar flux den- 
sity varied between 2 and 15 mJy while the RMS 
noise ranged from 0.1 to 0.9 mJy. 

The AIPS task JMFIT was used to fit ellip- 
tical Gaussians to the images in order to ob- 
tain astrometric positions and position errors, as 
listed in Table 1. The position errors listed for 
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each epoch have three components, which we es- 
timate and add in quadrature: (a) random er- 
rors in the position of the in-beam calibrator, es- 
timated as the beam full width at half-maximum 
(FWHM)/(2 x SNR), (b) random errors in the po- 
sition of the pulsar, estimated as the fitted Gaus- 
sian FWHM/(2 x SNR), and (c) systematic errors 
in the pulsar position, contributed by the residual 
ionosphere between the in-beam calibrator and the 
target, as well as other possibly unknown errors. 
The deviation of the observed pulsar image from 
that of a point source convolved with the dirty 
beam (or point spread function) provides a mea- 
sure of the residual systematic errors. We take 
the deconvolved size of the pulsar (0^, which is 
zero when the image is truly point-like) and use 
the quantity 6 d / y/ (N ant - l)T ohs /T iono as an esti- 
mate of this last error component. This quantity is 
not sensitive to the choice of Tiono, the ionospheric 
coherence time, which is empirically determined 
as ~ 0.15 hour. The number of antennas N ant is 
usually 10 (or 9 when we lost a VLBA station due 
to inclement weather) while the observation time 
T obs ranged from 1 to 3 hours. The small separa- 
tion (12') of the pulsar and the in-beam calibrator 
implies that any systematic errors in the correla- 
tor model (Ma et al. 1998) are negligible, and 
any residual systematic errors in the astromctry 
are likely to be small (< 0.1 mas). 

2.3. Likelihood Analysis for Astrometric 
Parameters 

To derive the best fit values and confidence in- 
tervals for the parallax and proper motion of PSR 
B0919+06, we calculate a likelihood function for 
the parameters that describe the apparent pulsar 
position as a function of time. These parame- 
ters include (Aao,A(5o), the offset of the pulsar 
from an arbitrary assumed position at the refer- 
ence epoch (1997.0); (^i a ,fis), the proper motions 
in right ascension and declination; and w, the an- 
nual trigonometric parallax. 

The likelihood function is calculated as follows: 
iv data ^ 

£ = II 9 ^ Ul ~ "OAOff - . (!) 

i=l 

where 

<Sj, Si = di,Si (A«o, ASq, us, tt; U) 



are model estimates and cm,8i are the observed 
values of the pulsar coordinates at epoch tj (offsets 
from the assumed position for 1997.0, as listed in 
Table 1), A^ata is the number of data points (10 in 
this case), and g is a normalized Gaussian function 
with zero mean and unit variance. 

To obtain the marginal probability distribution 
of each parameter e 9 = (Aa , A5o, fJ. a , Us, tt), 
we calculate the normalized integral of the likeli- 
hood function over all other parameters: 

UW - Jd9C(9) ■ (2) 

From the marginal distributions for each parame- 
ter, we obtain the best fit (median) values, as well 
as 68% (la) confidence intervals. These are sum- 
marized in Table 2, along with derived estimates 
for the distance D and transverse velocity V PJ _. 

Figure 2 shows the residual position offsets af- 
ter subtracting the best fit proper motion, with 
the best fit parallax sinusoids overplotted. The 
marginal probability distribution of the parallax is 
plotted as an inset in Figure 2. The parallax, tt = 
0.83±0.13 mas, differs significantly from Fomalont 
ct al. (1999), who report tt = 0.31 ± 0.14 mas us- 
ing a subset of the current data set (4 epochs from 
1994-1998). It is now apparent that the position 
errors were underestimated, leading to a parallax 
estimate 4a displaced from the current result. To 
test the robustness of the fit reported here, the 
likelihood analysis was repeated after eliminating 
data points singly (1994, 2000) and in pairs (1994 
and 2000, both points from October 1998). In all 
cases, the best fit parameters remained consistent 
(though with greater uncertainty) with the values 
for the fit including all data. For example, the 68% 
confidence interval for the parallax ranged from 
0.81 ± 0.15 mas to 0.87 ± 0.15 mas, compared to 
0.83 ± 0.13 mas when using all the data. 

We also note that the likelihood of zero par- 
allax is ^8 orders of magnitude lower than the 
likelihood of tt = 0.83 mas, and the null result 
is ruled out at over the 6a level. This result im- 
plies that even if B09 19+06 were twice as far away 
(<~ 2.4 kpc), the parallax signature of 0.4 mas 
could have been detected at 3a. With a suitable 
in-beam calibrator and the addition of more sensi- 
tive antennas (Arecibo, GBT) to the VLBA, a par- 
allax signature should be detectable to 10 kpc at 
5a with <~ 25 observations. Beyond this, the vari- 
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Table 1 

astrometry measurements for psr b09 19+06 



Pulsar Position Offsets 
Observation Date East-West (mas) North South (mas) 



1994.857 . . . 


-45.89 ± 0.64 


-263.20 ± 1.41 


1995.728 . . . 


-30.69 ± 0.33 


-189.60 ± 0.59 


1996.252 ... 


-22.20 ± 0.38 


-143.80 ± 0.76 


1998.233 ... 


13.86 ± 0.31 


27.70 ± 0.73 


1998.244 . . . 


14.60 ± 0.32 


28.30 ± 0.77 


1998.778 . . . 


26.52 ± 0.89 


75.40 ± 1.02 


1998.797 . . . 


27.27 ± 0.97 


77.20 ± 1.35 


1999.917... 


46.34 ± 0.20 


173.20 ± 0.52 


1999.947 . . . 


46.93 ± 0.17 


175.80 ± 0.38 


2000.307 . . . 


52.30 ± 0.49 


207.30 ± 1.03 



a Measured from 09 h 22 m 14?0011, +06°38'22'.'7180, referenced 
to in-bcam position 09 h 23 m 03?8991, +06°38'58'.'9980 (J2000.0). 



Table 2 

astrometric parameters for psr b09 19+06 



Parameter 


Median value 


la Error 


^ a (mas yr" 1 ) 


18.35 


0.06 


Us (mas yr _1 ) 


86.56 


0.12 


7r (mas) 


0.83 


0.13 


D (kpc) 


1.21 


0.19 


V p± (km s- 1 ) 


505 


80 
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able source structure of calibrators at the 0.1 mas 
level (Fey, Clegg & Fomalont 1996) places a basic 
limit on current astrometric techniques. 
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Fig. 2. — Residual position offsets in right ascen- 
sion and declination after subtracting the proper 
motion; sinusoids corresponding to the best fit 
parallax of 0.83 mas are overplotted. Inset: The 
marginal probability distribution for the parallax, 
log 10 / 7r (7r) plotted against n (mas). A null result 
is excluded by ~8 orders of magnitude, or over 6a. 

2.4. Transverse Velocity of PSR B0919+06 

From our measurement of the total proper mo- 
tion (88.48 ±0.13 mas yr" 1 ) and parallax (0.83 ± 
0.13 mas) of B0919+06, we can estimate its trans- 
verse velocity to be 505 ± 80 km s _1 . This is con- 
sistent with the mean pulsar population velocity 
of <~ 450 km s" 1 estimated by Lyne & Lorimcr 
(1994), as well as the two-component model with 
characteristic speeds of 175 and 700 km s _1 ob- 
tained by Cordes & Chcrnoff (1998). 

3. Interstellar Scintillation: Data and 
Analysis 

Electron density fluctuations along the line of 
sight through the interstellar matter (ISM) scatter 
radio signals. For pulsars, this phenomenon has 
been studied extensively through dynamic spec- 
tra (intensity variations as a function of time and 
frequency), whose characteristics allow inferences 



about the distribution of scattering material along 
the line of sight. Much of the information in dy- 
namic spectra can be condensed into two param- 
eters derived from the autocorrelation function 
(ACF): the scintillation timescale Aid (defined as 
the 1/e width of the ACF intercept on the time lag 
axis) and the decorrelation bandwidth A^d (the 
half-width at half maximum of the ACF intercept 
on the frequency lag axis). 

For PSR B0919+06, measurements of these pa- 
rameters exist over a 20 year timespan, from 1980 
to 2000 (Cordes, Weisberg & Boriakoff 1985; 
Cordes & Wolszczan 1986; Bhat, Rao & Gupta 
1999, and unpublished data from Arecibo Obser- 
vatory in 2000). Cordes & Wolszczan (1986) ob- 
served this pulsar in 1984-85, and report inter- 
ference fringes due to multipath propagation ef- 
fects in some of their dynamic spectra, signifying 
that the ISM had caused multiple imaging of the 
pulsar, while other epochs show random structure 
consistent with single-image diffractive scattering. 
Bhat, Rao & Gupta (1999) observed this pulsar 
during 1994, and report gradual and systematic 
variations in the dynamic spectra, but no fringing 
events. The derived scintillation parameters for 
the entire set of data are plotted in Figure 3 (top 
and middle panels): as expected, these parameters 
vary significantly with time. 

In the analysis that follows, we do not ad- 
dress either the statistical errors in quantifying 
Az/d, Atd due to the finite number of scintles ob- 
tained in each observation or the measurement un- 
certainties due to the finite number of indepen- 
dent experiments. These errors are adequately 
represented by the scatter in the data, since we 
treat each observation independently. Refractive 
effects may also modify Ai^d and Atd, particularly 
in the epochs where multipath effects are impor- 
tant. These effects are difficult to quantify, but we 
account for their presence in the interpretation. 

3.1. Scintillation and Scattering Formal- 
ism 

For the analysis, we adopt the theoretical 
framework of Cordes & Rickett (1998, hereafter 
CR98), who provide a general treatment of diffrac- 
tive ISS without assuming any specific scattering 
geometry (like a thin screen or uniform medium). 
This generality neccessarily leads to some alge- 
braic complexity, which proves worthwhile in the 
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Fig. 3. — Multifrequency scattering parameters 
for PSR B0919+06: (Top) decorrelation band- 
widths A^d; (Middle) scintillation timescales Aid! 
and (Bottom) scintillation speed Viss,5/3, u ob- 
tained over 20 years (see text for references). 
The pulsar speed derived from VLBI parallax and 
proper motion is plotted with a dashed line. 

analysis. 

The scintillation speed Viss is defined as the 
speed of the ISS diffraction pattern with respect 
to the observer, which includes contributions from 
the velocities of the earth and the medium as well 
as the pulsar. Given values for D (in kpc) and ob- 
serving frequency v (in GHz), and measurements 
of A^d (MHz) and Aid (sec), we obtain Vass in 
km s- 1 (CR98, eq. 13): 



VlSS, 5/3, u — ^ISS,5/3,u 



vAt d : 



(3) 



where A 



ISS, 5/3, u 



= 2.53 x 10 4 km s"\ and the 



subscripts 5/3 and u denote a Kolmogorov scat- 
tering medium that is statistically uniform. The 
estimated scintillation speed is plotted for the data 
set in the bottom panel of Figure 3. 

To estimate the transverse speed of the pulsar 
V P j_ from Viss 7 we assume that the speeds of the 
observer and the medium are small compared to 
the pulsar's; then (CR98, eq. 23): 



V p± = VF C VFd,PmVisS,5/3,u, 



(4) 



relate the scintillation pattern speed to the pulsar 
speed, and rescale the velocity from the uniform 
Kolmogorov medium case to the actual medium. 
Wq depends on the wavenumber spectrum of the 
medium as well as its distribution (CR98, eq. 18), 
and its contribution is small (maximum range of 
0.9-1.3) unless one assumes solely a thin screen for 
the distribution of the entire scattering medium. 
For the remaining analysis, we assume this quan- 
tity is unity. 

W^d.pm is a weighting factor determined by the 
distribution of the scattering material along the 
line of sight, which is proportional to the coef- 
ficient of the electron-density wavenumber spec- 
trum C^(s) (CR98, Appendix A), where s is the 
distance measured from the pulsar towards the ob- 
server, and varies from to D. An exact expres- 
sion for Wd,pm is derived for a square law phase 
structure function in CR98 (eq. 25): 



Wd.pmOD) 



2f D ds( S /D)(l-s/D)CUs) 



1 1/2 



J Q D ds(l-s/DyC^s) 



(5) 

This expression is only approximately valid for a 
true Kolmogorov medium, but the approximation 
is sufficient for our analysis. 

3.2. Hybrid Analysis of VLB and ISS Data 

There are two independent measurements of the 
transverse velocity of the pulsar, one from the ob- 
servable VLB proper motion fi, and the other from 
the measured scintillation parameters Af d , Ai d at 
an observation frequency v. 



v P± 



(6) 



lV C IVD,PMil SS ,5/3,u*. (7) 



where Wc and Wd,pm are weighting factors which 



Wd,pm depends on the (unknown) distribution of 
scattering material, as expressed in equation (5), 
while D is usually known only roughly from dis- 
persion measure-distance models (TC93). 

In general, if a parallax distance is unavailable 
for a pulsar, the different dependencies of V PJ _ 
on D can be used to iterate the two equations 
(6,7) for a given distribution of scattering mate- 
rial (uniform, exponential, TC93), and thus solve 
for the distance as well as the transverse velocity. 
This procedure is especially useful for pulsars well 
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above the Galactic z scale height, where dispersion 
measure-based models can provide only lower lim- 
its on the distance, and parallaxes are not likely to 
be obtained in the forseeable future (as discussed 
in §2.3). 

Alternatively, it is possible to refine electron 
density models by postulating the existence of ex- 
tra scattering material along the line of sight, in 
the form of a clump or a screen which contributes 
increments in dispersion measure ADM and scat- 
tering measure ASM at a location D s . Then 
a set of constraints on the acceptable values of 
the screen parameters (ASM, D s ) as well as the 
(poorly known) distance D can be derived from 
this hybrid analysis. 

3.3. Hybrid Analysis using the Parallax 
Distance for B0919+06 

For PSR B0919+06, the measured parallax dis- 
tance allows stronger constraints on the scatter- 
ing geometry. The observation of fringing events 
for this pulsar, as well as the wide variability in 
its scintillation parameters (Cordes & Wolszczan 
1986; Bhat, Rao & Gupta 1999), suggest that 
this pulsar is viewed through a localized screen or 
clump of variable scattering strength in addition 
to a distributed medium. Using only the TC93 
electron density model results in an overestimate 
of the distance to this pulsar, and fails to account 
for either the observed values of the scattering pa- 
rameters or their variation with time (Figure 3). 

We note that large changes in scintillation pa- 
rameters can be produced with relatively small 
changes in the average electron density at each 
epoch. From the definitions of dispersion measure 
and scattering measure (DM = J® dsn e {s); SM 
= ds C 2 (s)), we have: 

ADM = n e , s As, (8) 
ASM = Cl s As, (9) 

Cl s <x F s nls, (10) 

where F is a dimensionless "fluctuation parame- 
ter" as defined in Taylor & Cordes (1993), n e is 
the electron density in the medium, the subscript s 
designates the values for a screen of thickness As, 
and ADM, ASM are epoch-dependent values. In 
order to fit the observations, we postulate the ex- 
istence of a thin patchy screen of scattering mate- 
rial along the line of sight in addition to the TC93 



scattering material: 



C 2 n (s) = Cl TC93 (s) + ASMS(s ~ D s ), (11) 



where the thin screen contributes an additional 
scattering measure ASM = C 2 s As at a loca- 
tion s — D s , represented by the delta function 
S(s — D s ). Recasting equations (8-10) in terms of 
the differential DM contributed by the screen, and 
putting in the appropriate scaling factors, 



ASM 



(l/3)(2 7 r)- 1 /3 



K u F t 



(ADM) 2 
As 



(12) 



where the scale factor K u = 10.2 x 10~ 3 m~ 20 / 3 cm 6 
provides the appropriate unit conversion for ASM 
in kpc m~ 20 / 3 , ADM in pc cm -3 and As in par- 
sees. 

We can rewrite the weighting function VFd.pm 
from equation (5) in terms of the TC93 model and 
this additional scattering screen: 

Wt,,pm(D) EE 

2f Q D ds (s/D)(l - s/D)Cl TC9:i (s) + (D 3 /D)(l - D./D) ASM 



£ ds (1 - s/D)2Cl TC93 (s) + (1 - D./D)' ASM 

We integrate the TC93 model out to the paral- 
lax distance derived for this pulsar, which gives 
Wd,pm as a function of two screen parameters in 
equation (13), i.e. the relative scattering strength 
ASM/SM, and location D s /D. 

Additionally, from equations (6,7), we solve for 
the values of Wd,pm required to match the scintil- 
lation and interferometric velocity estimates (with 
Wq set to unity): 



(13) 

1/2 



4lSS,5/3,u 



D 
Aim 



1/2 



(14) 



Thus Wd,pm is known as a function of ASM/SM, 
D 8 /D (eq. 13), and the values of Wd,pm required 
for the two estimates of V p± (eqs. 6,7) to agree 
at each observation are specified by equation (14). 
This information is summarized in Figure 4, where 
the required values of Wd.pm are plotted as con- 
tours against ASM/SM and D s /D on the surface 
defined by equation (13). The range of values re- 
quired for the scattering strength and location of 
the screen (in addition to the TC93 model) can 
thus be determined. 
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Fig. 4. — Wd,pm expressed as a function of 
scattering screen strength ASM and distance D s 
from the pulsar towards the observer. The con- 
tours represent the value of Wd,pm required at 
each epoch shown in Figure 3 for the scintil- 
lation speed to match the actual V p , . Panels 
show slices through the surface on which contours 
are drawn, taken at the locations of the arrows. 
Top panel: ASM/SM = 0.1,10,1000 (thin to 
thick line); side panel: D s /D = 0.025, 0.050, 0.975 
(thin to thick line). Assuming a single screen, 
its location D s /D must be consistent with all of 
the contours (barring strong refractive scintilla- 
tion effects) , while the range of contours along the 
ASM/SM axis gives the required variation in the 
scattering strength of the screen. 

3.4. Results from the Hybrid Analysis 

A reasonable assumption is that a turbulent 
screen may fluctuate in scattering strength on 
short timescales, while its location varies slowly (if 
at all). In that case, the screen must be located 
at a distance consistent with all of the observed 
contours in Figure 4, while the spread in the con- 
tours along the ASM/SM axis gives a measure of 
the range over which the screen scattering strength 
has varied over the 20 years spanned by observa- 
tions. From the figure, consistency with the inner- 
most contours requires D s /D > 0.9. However, the 
multiple imaging and fringing events observed dur- 
ing 1986 indicate the occurrence of strong refrac- 
tive scintillation effects. The current analysis of 



diffractive scintillation probably overestimates the 
scattering required during that epoch. Thus the 
innermost set of contours in Figure 4 may not be 
useful in locating the screen, and the constraints 
on the screen can be relaxed to D s /D > 0.8, with 
ASM/SM ranging from 0.1 to 100. 

For a pulsar distance of 1.2 kpc, this result sug- 
gests that the scattering screen is within ^240 par- 
sees of the Earth, implicating structure within the 
local ISM. Assuming a screen thickness As = 10 
parsecs, equation (12) gives an approximate range 
for ADM and the fluctuation parameter F s : 

10~ 5 < F S (ADM) 2 < 10~ 2 (pc cm" 3 ) 2 . (15) 

Assuming, for example, that F s = 10 (consistent 
with values in TC93, Table 2), the required range 
in ADM is 10~ 3 to 3 x 10~ 2 pc cm" 3 . For the as- 
sumed screen thickness of 10 parsecs, this requires 
rt e , s between 10 -4 and 0.003 cm -3 . These values 
are consistent with the DM variability observed for 
several nearby pulsars (ADMrms ~ few x 10~ 3 
pc cm" 3 ), including B0919+06 itself, which had a 
DM variation of 4 x 10" 3 pc cm" 3 in 1989-1991 
(Phillips and Wolszczan 1992). Of course, in real- 
ity, both F s and the screen thickness As may vary 
along with changes in electron density, but these 
values show that the derived constraints on screen 
location and ASM/SM are well within the range 
of physical possibilities. Additionally, simultane- 
ous measurements of DM variability and changes 
in ISS parameters can be used to estimate As, 
which is related to the outer scale size of ISM fluc- 
tuations. 

3.5. Interpreting the Screen in terms of 
the Local ISM 

Combining the analysis for PSR B0919+06 
with parallaxes and scattering measurements for 
other pulsars in the third Galactic quadrant pro- 
duces a coherent picture of the local ISM in this 
direction. Measurements of Galactic coordinates 
(i,b), distance, and DM for third-quadrant pul- 
sars are summarized in Table 3, along with de- 
rived values of n e , angular separation A8 and 
Galactic z height. Toscano et al. (1999a) review 
the local ISM, and their Figure 2 provides use- 
ful background for our discussion. The preceding 
analysis, and a recent (upward) revision in the 
distance to PSR B0950+08 (Brisken et al. 2000), 
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allows us to improve upon their discussion for the 
third Galactic quadrant. 

It has long been known that there is an elon- 
gated cavity in the local neutral ISM, surrounding 
the Sun and extending several hundred parsecs to- 
wards Galactic longitude I = 240°. This feature 
appears to be correlated with the region of low 
reddening in the Gould Belt, between t = 210° 
and 255° (Lucke 1978), and the much- reduced 
H I absorption at both high and low latitudes 
in this longitude range (Paresce 1984). The 
Local Bubble fills part of this cavity out to a 
distance of 0.1-0.2 kpc (e.g., Cox & Reynolds 
1987), though the Bubble boundary is not well 
defined in this direction. Beyond the Local Bub- 
ble, Heiles (1998) has proposed a superbubble, 
GSH 238+00+09, along £ = 238°. This extends 
from 0.2 to 1.3 kpc from the Sun, and fills another 
part of the elongated cavity. 

It is expected that the region between these two 
structures has clumps of dense and partially ion- 
ized gas, possibly left over from the creation of 
either (or both) of the bubbles, or from their in- 
teraction with each other and their environments. 
For example, Dupin & Gry (1998) report the ex- 
istence of two dense clouds of H II between 40 and 
90 pc from the sun, along the line of sight towards 
(3 Canis Majoris {£, b = 226°, -14°), and speculate 
that these originate in the formation of the Local 
Bubble. Heiles (1998) postulates "filaments" of 
warm ionized material on the surface of superbub- 
bles. The scattering properties of various pulsars 
toward the third Galactic quadrant (Table 3) are 
consistent with the existence of such a turbulent 
interface region. 

The most striking feature of Table 3 is the 
difference in n e between the lines of sight to 
PSRs B0950+08 and B0919+06. Although sep- 
arated by only 7? 8 (35 pc at a distance of 250 pc), 
PSR B0919+06 exhibits strong scattering and has 
undergone episodes of multiple imaging and re- 
fractive fringing, while PSR B0950+08 has one of 
the lowest scattering measures known (Phillips & 
Clegg 1992). There are two possible causes for the 
difference in the properties of these two pulsars. 

One possibility is that PSR B0950+08 may be 
closer than the turbulent interface region. This 
model requires the distance to the interface region 
to be larger than 280+25 pc, the (revised) distance 
to PSR B0950+08, or equivalently, D s /D < 0.8. 



We concluded above that D s /D > 0.8, but given 
the uncertainties in incorporating the effects of 
strong refractive scintillation into the estimate for 
D s /D, we can not rule out this possibility. Rather, 
we take D s /D ~ 0.75-0.8 as a reasonable range for 
the screen distance. 

An alternative possibility is that PSR B0950+08 
is within or beyond the turbulent interface re- 
gion, but the turbulence is concentrated into 
clumps. In order to cover PSR B0919+06 but 
not PSR B0950+08, the typical (transverse) scale 
size of a clump would have to be less than roughly 
30 pc. This suggestion is consistent with the 10 
pc (radial) scale size of the scattering screen that 
was assumed in our analysis. 

The lines of sight to pulsars B0823+26 and 
J1024— 0719 also provide useful information, 
though the derived constraints are weaker. Both 
pulsars have mean electron densities along the line 
of sight comparable to or larger than B0919+06. 
Both exhibit stronger scattering than B0950+08, 
particularly B0823+26 which has undergone 
episodes of strong refractive fringing and multi- 
ple imaging (e.g., Clegg, Fiedler & Cordes 1993). 
J1024-0719 is closer to the Sun than B0950+08, 
perhaps indicating that B0950+08 is within or 
beyond the turbulent interface region. However, 
both B0823+26 and J1024-0719 are near the 
edges of the H I cavity. Whether their scattering 
properties and mean electron densities are due to 
the "edge" of the H I cavity or to the interface be- 
tween the Local Bubble and the GSH 238+00+09 
superbubble is not clear. 

We conclude that interstellar scintillometry of 
PSR B0950+08 and PSR B0919+06 indicates a 
turbulent interface region between the Local Bub- 
ble and the GSH 238+00+09 superbubble. This 
turbulent interface region begins 250-300 pc from 
the Sun in the direction (£,b) w (225,40) and 
possibly contains clumps of scale size < 30 pc. 
This result is not inconsistent with the model 
of Bhat, Gupta & Rao (1998). Given the 
height of PSR B0950+08 above the Galactic plane 
(~ 200 pc), it is possible that the interface re- 
gion, and therefore the Local Bubble and the 
GSH 238+00+09 superbubble, extend to higher 
Galactic latitudes. 
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Table 3 

Pulsars in the Third Galactic Quadrant 



PSR 


B0919+06 


B0950+08 


J1024-0719 


B0823+26 


en 


225.42 


228.91 


251.70 


196.96 


b(°) 


+36.39 


+43.70 


+40.52 


+31.74 


DM (pc cm" 3 ) 


27.31 


2.97 


6.50 


19.48 


D (kpc) 


1.2 1 


0.28 2 


< 0.23 3 


0.38 4 


z (kpc) 


0.72 


0.19 


< 0.15 


0.20 


n e (cm -3 ) 


0.023 


0.011 


> 0.029 


0.055 


A6 (°) 





7.8 


21.4 


23.7 


Alj_ = D S A6 (pc) 





35 


93 


104 



References. — Distances from (1) this work; (2) Brisken et al. 2000; (3) 
Toscano et al. 1999b (upper limit derived from /i, P, and P); and (4) Gwinn 
et al. 1986. 



4. Conclusions 

Parallax distances to pulsars have the prospect 
of resolving several outstanding questions, both 
about the pulsar population and about the in- 
tervening ISM. In this work, we have presented 
VLBA astrometry on PSR B0919+06, using 
phase-referenced observations with an in-beam 
calibrator only 12' from the pulsar. The use of 
the in-beam source reduces ionospheric effects and 
other astrometric errors to the 0.1 mas level and 
allows sub-milliarcsecond astrometry. With ten 
observations (seven distinct epochs) distributed 
over seven years, we derive a proper motion 
fi a = 18.35 + 0.06 mas yr" 1 , ^ s = 86.56 + 0.12 mas 
yr" 1 , and measure a parallax 7r = 0.83 + 0.13 mas 
(68% confidence intervals). This result implies a 
transverse speed of 505 ± 80 km s" 1 at a distance 
D = 1.21 ± 0.19 kpc, making PSR B0919+06 one 
of the farthest objects for which a trigonometric 
parallax has been obtained. The use of an in- 
beam calibrator with the VLBA at 1.5 GHz can 
provide trigonometric parallaxes to twice this dis- 
tance, or even further (to 10 kpc) with the use of 
more sensitive antennas and higher frequencies. 

For the measured distance of 1.2 kpc and a 
DM of 27.31 pc cm" 3 (Phillips and Wolszczan 



1992), the derived mean electron density towards 
B0919+06 is 0.023 cm" 3 . This line of sight is along 
the same Galactic longitude as the superbubblc 
proposed by Heiles (1998), as well as clouds of 
ionized gas in the local neighbourhood (Dupin & 
Gry 1998), though at a higher Galactic latitude. 
We use scintillation data spanning 20 years in a 
hybrid analysis, equating the derived ISS and in- 
tcrfcrometric speeds of the pulsar in order to con- 
strain the distribution of scattering material along 
the LOS. We find that excess scattering mate- 
rial is required compared to the standard TC93 
model, and assuming that this excess material is 
confined to a thin screen or clump, we constrain 
its location to within ~ 250 pc of the Sun. Com- 
parison with the neighboring lines of sight to pul- 
sars B0950+08, J1024-0719, and B0823+06 in- 
dicates that this result is consistent with a tur- 
bulent interface region between the Local Bubble 
and the GSH 238+00+09 superbubble, which pos- 
sibly contains clumps of scale size < 30 pc. We 
note that using an updated model of the Galactic 
electron density distribution (Cordes & Lazio, in 
preparation) will not materially alter these con- 
clusions. 
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